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Abstract 

We investigated the Dirac electrons transmission through a potential barrier in the presence of 
circularly polarized light. An anomalous photon-assisted enhanced transmission is predicted and 
c/3 \ explained in a comparison with the well-known Klein paradox. It is demonstrated that the perfect 

transmission for nearly-head-on collision in an infinite graphene is suppressed in gapped dressed 
states of electrons, which is further accompanied by shift of peaks as a function of the incident 
angle away from the head-on collision. In addition, the perfect transmission in the absence of 
potential barrier is partially suppressed by a photon-induced gap in illuminated graphene. After 
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transmission with no potential barrier becomes completely suppressed and the energy range for 
the photon-assisted perfect transmission is reduced at the same time. 



I. INTRODUCTION 

Research activity on the properties of graphene has been growing rapidly ever since its 
experimental discovery and demonstration of its unusual properties arising from its energy 
band structure [l|, ^ . The novel properties of graphene may be attributed to its massless 
Dirac fermions at the Fermi energy J3|]. An interesting consequence of the Dirac electron is the 
Klein paradox j^ in which an electron in graphene undergoes unimpeded tunneling through 
potential barriers of arbitrary height and thickness. This property of Dirac electrons is due 
to their linear energy dispersion relation or helicity. Electrons are said to be chiral if their 
wave functions are eigenstates of the chirality operator h = a ■ p/(2p) where a = {a^, Cy} 
is the Pauli vector consisting of Pauli matrices and p = {p^^, Py} is the electron momentum 
in graphene layers. Electrons in graphene near the K points (around the corners of the 
hexagonal Brillouin zone) are chiral due to the fact that the chirality operator is proportional 
to the Dirac Hamiltonian which automatically makes chirality a good quantum number. 

The fact that the Klein paradox is also obtained in bilayer graphene makes this effect even 
more sophisticated. This leads us to realize that the Klein paradox is not simply due to linear 
electron dispersion but may be observed for both massless and massive quasiparticles j^ . In 
this paper, we consider a sharp p — n — p junction or potential barrier profile. This type 
of potential can be constructed by underlying n,etal contact ot insulattng stnp fl and was 
employed to demonstrate unimpeded transmission p|] . 

Both the tight-binding model and k ■ p approximations for infinite graphene sheet ac- 
curately show that electrons and holes have linear energy dispersions e{k) = shvp k = 



s^t^i^y/AjJ^r^ with no gap, where s is the electron- hole parity with s = 1 for electrons and 
s = — 1 for holes. For the potential barrier described above, there is a translational sym- 
metry in the y direction parallel to the boundaries of the potential so that ky is conserved. 
In contrast, the longitudinal component k^ is modified by the potential so that when the 



particle has energy E^ we have k^i = \/{E — ViY /{fivpY ~ ^yi where Vi is the potential in 
the region i. 

The electron effective mass and its properties will play a crucial role in our analysis. In 
infinite intrinsic graphene, the Dirac electron is massless although the mass can be intro- 



duced 



6|] or implemented experimentally. For bilayer graphene, the particle effective mass 



exists in any possible approximation [7| , leading to the existence of evanescent terms in the 
wave function j^ although this does not violate chirality symmetry and the Klein paradox. 

The outline of the rest of this paper is as follows. In Sec. HIl we derive the eigenstates and 
transmission coefficient for Dirac electrons dressed with photons from a circularly polarized 
light. Section UTTj is devoted to a formalism in which the roughness of the boundaries for the 
potential barrier is included phenomenologically in calculating the transmission coefficient 
and numerical results are presented. A brief summary is given in Sec. IIV[ 

II. TRANSMISSION COEFFICIENT FOR DRESSED ELECTRON STATES 

It was shown recently [8|, |9|] that when Dirac electrons in a single graphene layer are 
interacting with an intense circularly polarized light, electron states will be dressed by 
photons. The main idea of the present study is to investigate the transmission properties of 
such dressed electrons for the case of single layer graphene. We go beyond the approximations 
used in Ref. J8| by retaining the results up to the order of 0{A^) so that we are able to 
investigate the difference between the dressed states and massive Dirac electrons described 
below by the Hamiltonian in ( 124|) . Here, A is a quantity measuring the induced gap between 
the valence and conduction bands of dressed electrons. 

We begin with the electron-photon interaction Hamiltonian 

n = VFa- (p- e Acirc) , (1) 

where vp is the Fermi velocity and the vector potential for circularly polarized light of 
frequency uq can be expressed as 

Acirc = \ 77 (e+a + e_a^) = W [{a + a^)e^ + i{a - a^)ey] (2) 

V Co Wo v V ^eoi^ot/ 

in terms of photon creation and destruction operators a^ and a^, respectively. Here, V is the 
mode volume of an optical field. In order to study the complete electron-photon interacting 
system, we must add the field energy term tvjjQO^a to the Hamiltonian ([1]). As usual, we 
seek the wave function in the form of a plane wave \&(r) = e*''''"^(A;), which results in the 
following reduced Hamiltonian 



1-i = fvjjQ a^a + hvpcr ■ k 
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. The reduced Hamil- 



tonian ([3]) for infinite graphene can be understood as the sum of two parts, namely, the 
Dirac Hamihonian 



'Hoirac = hvp 0" ■ k = hvpicTxkx + CFyky) = hvF(a-k^ + cr+/i;_) 
and the Jaynes-Cummings Hamiltonian 



(4) 



Hj-c = hjjQ o)ci — 



w 
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(cr+a + cr„a^) 



(5) 



which corresponds to a two-level quantum optical system and, most importantly, can be 
solved analytically. Here, we have defined k± = k^ i: iky and A^ represents the number of 
radiation quanta (intensity) for the incident optical field. We only consider the situation 
such that the electron-photon interaction amplitude w is much less than either the photon 
or Dirac electron energy, i.e.. 



^2Nhe'^v1 



w 



2a huo <^ Jtwq , 



(6) 



where a = w/hwQ <^ 1. The two eigenstates of the Hamiltonian (^ could be obtained with 
an expansion over the basis of just two functions |l)Ar = | '[', A^) and |2)Ar = | |, A^+ 1) for 
each A^ value, that is, 



l^t.Af) = /iAf|l)Af + '^n\2)n , 
|*i,7v) = fJ'N\2)N - iyN\l)N , 

which corresponds to the following Hamiltonian in the chosen basis |l)Ar and |2) 
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(9) 



In this way, the transformation ([7]) and ([8]) becomes just a simple rotation in the Hilbert 
space. This expansion yields the eigenvalue equation and its solution below 



{e-Nnwo)[e-{N+l)nwo]-{^nuJoJ^^^] =0, (10) 



where e^- and ej, correspond to the lower + and upper — signs in the solution. For simplicity, 
we assume here the radiation is classically strong with A^ ^ 1. However, we note that we 
may set A^ + 1 ^ A^ only in the terms (9(a^) but not in the terms with NHuiq. Furthermore, 
it is a simple matter to obtain the expansion coefficients fiN and un as 

o 
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pn = sin 9c , 

aJ(N+l)/N a , , 

tan^c = ^ - • (12) 

l + v/l + «'(iV + l)/iV 2 ^ ' 

One can easily verify that the wave functions ([7]) and ([8]) are the eigenstates of the Jaynes- 

Cummings Hamiltonian (|5]) with energies (ITT]) , to any order of a. Here, for all the above 
derivations, we have employed the standard relations, i.e.. 



a}\ ti, A^) = ViVTTl ti, AT + 1) , a\ ti, A^) = ViV| ti, iV - 1) , 

a+|;,A^) = |t,iV), a+|t,iV) = 0, 

a_|t,iV) = |i,Ar), a_|i,A^) = 0. (13) 

We now look for the eigenstates of the Hamiltonian ([3]) as expansions over the set of Jaynes- 
Cummings Hamiltonian eigenfunctions (ITj) and (jH]). We will confine our attention to the 
field source with only three nearest photon occupation numbers, i.e., A^ = A'o — 1? ^o and 
A^'o + 1, leading to 

ATo + l 

\^{k)) ^ Y. (^M(fc) I *t,^) +C2,.(A:) I vl.;,.)) . (14) 



We know that the eigenfunctions ([7]), (j8]) corresponding to different numbers A^ are or- 
thogonal to each other. First acting the Hamiltonian ([3]) on f fl^ . and then, multiplying 
both sides of the expansion ( IT^ by (^i^at^I, (\&4^_7Vq_i|, (\&|^Ar,)+i| and (\&4^Ar(,|, this result in 
the following four equations 



Cl,No I ^0 ~ 



a^ Nq + I 



huo + hvp {^^,No\o' ■ k|$(A;)) = eCi^ 



No 



C2,No~i {No + — ^ _ ) hwo + hvF {^i,No~iW ■ k|$(A;)) = eC2,No~i , (15) 
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a^ No + 1 



hwo + hvp (*t,Afo+ik ■ k|$(/i;)) = eCi^No+i , 
hwo + hvF (^i,jVok ■ k|$(A;)) =eC2,No ■ (16) 



4 No 

For chosen number A'o, each pair of these equations describe two energy subbands that 
are separated by ^ a'^huo/2 a.t k = 0. Here, we include only four nearest energy subbands, 
corresponding to the dressed states with different photon occupation numbers A^ and electron 
states with subband indices t (lower energy) and J, (higher energy). 

Taking into account the following relations for the Dirac Hamiltonian, 



a ■ k |^t,A^o) = f^Nok+\ i, No) - i^Nok-\ t, ^^o + 1) , 
a ■ k |^|,^„) = -fiNok-\ t, No + l)- UNok+\ i, No) 



(17) 



and with the simplifications described above, we can explicitly write out the Dirac Hamil- 
tonian terms in equations (IT5l) and flT6|) . The presence of photon occupation numbers does 
not follow from the dressed states Hamiltonian and cannot be determined only from the 
ratio between the photon filed energy Nfiwo and the electron-photon interaction amplitude 
w. It should results from the model of the circularly polarized light source. If we consider 
electrons near the K point such that the kinetic energy hvp k ^ a'^hwQ/2, we can only retain 
only three photon occupation numbers A^o — 1? ^o and A'o + 1, which is consistent with the 
approximation described above. Consequently, we arrive at the system which consists of 
weakly coupled equations to determine the two nearest subbands, i.e.. 
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where /i = fiNo, v = un^, A ^ a^hwo/A, Ci = Ci^No, C2 = C2,No~i, C'a 
C4 = C2,7Vo- This leads to the following energy dispersions with /i ^ 1: 



Ci 



(18) 
Afo+i and 



eiAk) = (no + ^)hiuo-lnwo^l + V + ^(1 + i^')k' ± 2y^^k^ + l)(eP + r]) , (19) 
esAk) = (no + ^]hLUo + \ nuJo\ll + V + i{l + ^^)k^ T 2^{^^k^ + l)(eP + r^) , (20) 

where ^ = {2vp/ujqY and 77 = (2A/^(X'o)^ ~ a^/4. Another approximation which we may 
use here is to consider only small wave vectors k so that we get two independent pairs of 
subbands separated by energy hujQ. This gives a two-component spinor wave function where 
each component consists of two independent terms. If we consider the transmission of such 
states through a potential barrier whose height is equal to /Iojq, one of the terms in region 
"1" will exactly match the other in the region "2" (the potential region) so that a part of 
the wave function will be completely transmitted. This will lead to a substantial increase in 
the total transmission. On the other hand, if we consider a larger number of such subband 
pairs, this effect will not be significant since this complete transmission occurs only for one 
of the wave function terms. This is the reminiscent of the Klein paradox since it does not 
depend on the barrier width as long as the barrier height is exactly equal to fvjjQ. Revealing 
such an anomalous increase for the dressed states tunneling is an important discovery of the 
present paper. The corresponding eigenvalue equation for the case of two independent pairs 
of subbands yields, as expected 



{h\l e + A^-^e- No hujof} {ffvl e + A^-ie- {No + 1) hujof} = 



(21) 



The simplest possible approximation is to keep only two terms with the coefficients 
[Ci(/c), C2{k)] leading to the two nearest energy subbands separated by a gap 2 A ^ a'^huo/2 



at fc = due to electron-photon interaction. This approximation results in a simplified 
algebraic system determining the pair of coefficients Ci{k) and C2{k), i.e., 



(A^o ^0 — ^)Ci + hvp {kx + iky) C2 = eCi 
{No fkuo + A) C2 + hvp {kx - iky) Ci = e C2 



(22) 



The non-trivial so 
obtained in Ref. \t 



ution of these equations gives the energy dispersion which was previously 

1 as 



e{k) = NohuJo± y'A2 + h^v^ k^ = Nohjo + f3^A^ + h^v^ k^ . (23) 

If the electron-photon interaction is removed, then a = w/2 -^ 0, and the energy dispersion 
relations fl23|) demonstrate a non-interacting system consisting of a Dirac electron ep{k) = 
Phvplkl and photons NofiuiQ. All the other energy subbands are separated at least by 
fiiJo ^ w and could be neglected, which justifies the above two-subband approximation. In 
this notation, /3 = ±1 is the dressed conduct ion /valence band index corresponding to bare 
electron/hole bands for infinite graphene when A — t- 0. 



The system fl22]) is formally similar to the eigenvalue equations for the case of the effective- 
mass Dirac Hamiltonian 



'H = hvFa-k + V{x) 



1 
1 



A (73 



(24) 



where a^ is a Pauli matrix and V(x) is a one-dimensional potential. The electron dispersion 
and transmission pro per ties for both a single as well as multiple square potential barriers 
have been studied [y, [lO| for monolayer and bilayer graphene [ll| . It was also shown that 



a one-dimensional periodic array of potential barriers leads to multiple Dirac points 12 1. 
Several papers have introduced an effective mass term into the Dirac Hamiltonian for infinite 



graphene which may be justified based on different physical reasons 13|. For example, it 
has been shown [ij] that an energy bandgap in graphene can be created by boron nitride 
substrate resulting in a finite electron effective mass. However, we emphasize that the 
analogy between the Hamiltonian ( !24l) and that for irradiated graphene is not complete since 
that would correspond to A < 0. Although this difference does not result in any modification 
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of the energy dispersion term containing A^, it certainly modifies the corresponding wave 
function. 

The interaction between Dirac electrons in g raphene and a circularly polarized light has 



been considered in the classical limit in Ref. 



15| . In this limit, a gap in the Dirac cone opens 



up due to nonlinear effects. The dressed state wave function has the form 



^dr{k) = 



with Ci{k) 7^ C2{k) given by 



Ciik) 
(3 C2{k) e'^ 



(25) 



Ct{k) 
CHk) = 



2(1+72)^27/1+^ 



v/1 + 7^ + 7 



(26) 
(27) 



^2(l + 72)^27v/l + 7' 

corresponding to the energy subbands e{k) = Nq hujQ ± y^A^ + h?Vp k"^. In this notation, 
7 = A/{hvF k) and is the angle which k makes with the longitudinal x axis. Without an 
optical field, i.e. A = 0, we obtain Cf = C^ = l/-\/2. In the limit (A ^ A;), the coefficients 
exhibit peculiar symmetry with 



C7,{k) 
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+ 



37=^ 



V2 2V2 8^2 16^2 ' 

^i ^^^+ 37=^ 



(28) 



V2 2^2 8^2 16^2 ' 
We note that this expansion is not valid too close to the Dirac point and should not be used 
for arbitrary wave vector to calculate, for example, the polarization function. Additionally, 
one may verify that Ci{k) 7^ C2{k) for any chosen A in the range of validity. Consequently, 
the chiral symmetry is broken for electron dressed states 



h'i/drik) 



1 O" ■ p 

2 p 



PC2{k) 

Ci{k)e'^ 



(29) 



cm 

Clearly, it follows from ( 129|) that the non-chirality of the dressed electron states becomes 
significant if the electron-photon interaction (the leading 7 term) is increased. This affects 
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the electron tunneling and transport properties. We now turn to an investigation of the 
transmission of electron states through a potential barrier when graphene is irradiated with 
a circularly polarized light. 

For simplicity, we consider a square potential barrier of height Vq given by V(x) = 
Vq [0{x) — 6{x — Wq)] where Wq is the barrier width and 6{x) is the Heaviside step func- 
tion. Since the wave number kx is the same in region 1 (x < 0) and region 3 (x > Wq) 
and the current component is jx = $Va$, we only need the wave-function continuity at 
the potential boundaries for the system considered. From this continuity condition, the 
transmission probability T can be determined from T = |tp where t is the transmission 
coefficient or the amplitude of the wave propagating forward in region 3. Here, we only 
show an analytical expression for the transmission coefficient of the dressed states in the 
limit of 5 ^ Vo corresponding to kxi ^ kx2'- 



cos^0 3 cos^ {kx2Wo) cos"^ (p sm^ {kx2Wo) 2 



cos^ {kx2Wo) cos^ + sm'^{kx2Wo) ffvj, kl^ [sin^(A;^.2Wo) + cos'^{kx2WQ) cos^ 

(30) 

Here, 6 = tan~^{ky/kx2) — > 0, </> = ta,n~^{ky/kxi), the second term includes the effect 
of electron-photon interaction (oc A^) to the leading order. In addition, we only show 



those terms of the lowest order in e/Vo. There is another relevant study [16|], investigating 
the tunneling of Dirac electrons with a finite effective mass, a parabolic dispersion in the 
presence of a energy gap, and a certain chirality, through a potential region. In that study, 
the particle tunneling through a square potential barrier differs from both Dirac electrons 
and the dressed states of electrons under a circularly polarized light illumination. 

For nearly-head-on collision with ky <^ kxi ^ kx2 for high potential as well as for infinite 
graphene (A — )■ 0), transmission coefficient has the following simplified form 

T = 1 - sin^{kx2Wo) {e^ - 2/39(P + 0^) , (31) 

where we assume Vq ^ e, 6' ^ ^ 1 and /3 = ±1. 
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III. NUMERICAL RESULTS AND DISCUSSION 

In our numerical calculations, energies will be measured in units of {?)kpat/2) with the 
carbon-carbon distance a ~ 1.42 A and the hopping parameter t = 2.7a/3/2 eV. We 
measure the wave vector in units of the Fermi wave number kp and write its components as 
kx = cos (j) and ky = sin in terms of the angle of incidence (j). 

In Fig.lH we present the transmission for dressed electron states with arbitrary energy 
and angle of incidence. We clearly see that dressing ruins the Klein paradox in (a) for 
head-on collision with = 0. The resonant peaks are shifted for the other incoming angles 
in (b) and the effect is stronger for small incident angles. In (c) and (d), the transmission 
probability plots are given in terms of the longitudinal momentum k^i in front the barrier 
and kx2 in the barrier region. We find that the intensity and locations of the transmission 
peaks in (d) are distorted compared to infinite graphene in (c). The diagonal kxi = kx2 
corresponds to the absence of potential barrier and should yield a complete transmission 



for A = 0. However, the condition {^/{e — Vq)^ — A^ > hvpky) for dressed states must be 
satisfied, which makes the diagonal incomplete (missing diagonal for small kxi and kx2) due 
to the occurrence of an induced gap. 

Figure [2] displays the effect due to electron-photon interaction on the electron transmis- 
sion in terms of incoming particle energy e and angle of incidence (p. From the figure, we 
see the Klein paradox as well as resonant tunneling peaks in the transmission probability 
for regular infinite graphene with A = in (a). The dark "pockets" on both sides oi e = Vq 
in (a) demonstrate zero transmission for the case je — Vo| ^ e, which results in imaginary 
longitudinal momenta kx2 for most of incident angles and produces a completely attenuating 
wave function. When a small gap is opened in (b) for dressed states of electrons in graphene 
under the illumination by a circularly-polarized light, we observe a set of complete trans- 
mission branches, where a strong dependence on </> for lower branches is seen. However, this 
dependence is greatly suppressed when the dressed-state gap is increased in (c), leaving 
us a set of equally-distant branches due to photo-assisted electron tunneling. 

As mentioned above, by including more than the two nearest subbands, as shown in in 
Fig.Et^a), the electron dispersion, wave function and transmission amplitude will be modi- 
fied. According to the approximation adopted by Kibis |8[, the two subband pairs may be 
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considered as independent. Under the condition of equal transverse momentum ky for both 
terms of the particle wave function, the first term of the wave function in region "2" with a 
potential barrier is similar to the second term in the regions without potential. Therefore, 
the states exactly match across the potential boundary, which should definitely increase the 
total transmission amplitude. The other possibility is that incoming angle and momentum 
of the second term is totally independent of the first one and leads to resonant transmission 
regardless of the transmission amplitude of the first term. Based on our derived results, 
we find that the transmission should increase even for non-split energy subband pairs and 
corresponding wave function. Since the existence of certain photon occupation numbers A^ is 
determined by the laser source, we can consider only one pair of the dressed states subbands 
is occupied while the other subband pairs are unpopulated. As illustrated in Fig.|3I^b), in the 
potential region "2" , a particle may populate another subband corresponding to a different 
number of photons instead of changing its longitudinal momentum for the barrier heights 
exactly equal to the multiple of /zwq. The opposite transition will occur at the boundary 
between regions "2" and "3". This will result in unimpeded tunneling T = 1 independent 
of the barrier width, as seen from Fig.[3]^c), which is expected to be a major contribution to 
the current. 

We now investigate the effect of disorder on the transmission probability through a poten- 
tial barrier in graphene. This can appear as short-range disorder, inter-valley scattering and 
trigonal distortion. In a single layer graphene, disorder also induces a metal-insulator tran- 



sition by creating a dynamical gap [17[ • Consequently, this effect can modify the gap created 
by the electron-photon interaction. In bilayer grpahene, disorder directly leads to energy 
dispersion with a gap as well as modify the energy dispersion close to the band edges. Addi- 



tionally, disorder may also lead to localized states inside a gap in bilayer graphene 18|, Il9 |. 



In this paper, we introduce disorder phenomenologically through the non-conservation of 
the transverse electron momentum ky. As mentioned above, ky is conserved for both Dirac 
electrons in infinite graphene and dressed states by photons. Introducing the quantity F 



as a measure for the disorder 
yielding 



20 



21| . we model its stochastic distribution as a Lorentzian, 
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U^s{k., ky) - - j^ dqy ^^-^-il|L , (32) 

where to{kx, Qy) denotes the transmission coefficient in the absence of any imperfections. 
As long as the disorder is weak with F ^ 1, different distributions, which give 5- function 
in the hmit of F — > 0, will result in almost equal transmission coefficients. For nearly- 
head-on collision, ky <C kx2, the transmission coefficient may be obtained analytically using 
a Gaussian distribution. The imperfect boundary of the potential region can also be the 
result of some stochasticity of kx2 to make the effect stronger. We neglect this effect since 
our goal is to investigate the role played by disorder using a simple approximation. 

Our numerical results showing the effect due to disorder are presented in Fig. HI First 
we test our numerical results in (a) by applying the Lorentzian transformation in (1321) to 
the complete transmission with T = 1, corresponding to the Klein paradox with the head- 
on collision for infinite graphene. Analytical integration clearly results in the transmission 
amplitude equal to unity. The transformed distribution demonstrates the precision of our 
numerical procedure. By comparing Fig.Hl^b) with Fig.[T]^d), we see a complete suppression 
of the perfect transmission with T = 1 along the diagonal kxi = kx2 by disorder along the 
boundary as well as a reduced range of kx2 for photon-assisted perfect transmission. 

IV. CONCLUDING REMARKS 

According to recently published results |8| , the interaction between Dirac particles in 
graphene and circularly polarized light leads to the formation of quantum electron dressed 
states. These states are appreciably different from conventional Dirac electrons in ordinary 
infinite graphene. From Ref. J8[, laser power 100 mW leads to a gap on the order of A ^ 
100 meV, which is required to make the effect significant for infrared light frequencies and 
room temperature. This enables possible experimental demonstrations of the described 
effects. We have shown that electron-photon interaction gives rise to states of broken chiral 
symmetry. The non-symmetrical properties of the states become more significant when 
the electron-photon interaction is increased. In addition, there are no dressed states with 
chirality symmetry. In general, incoming electrons or holes passing unimpeded through a 
square potential barrier require chiral symmetry. Under the illumination from a circularly- 
polarized light, we can control the degree of partially-broken chiral symmetry in dressed 

13 



states or the degree of partially-perfect transmission through a potential barrier in a graphene 
layer. 

In the simplest approximation when only the two nearest subbands are included, the 
model is formally similar to the so-called a^ Hamiltonian used to describe the particles in 
a single layer graphene with parabolic energy dispersion, giving non-zero electron effective 
mass. We discussed the similarities as well as the differences affecting the wave functions 
but not the energy dispersion. By including more the next-nearest subbands, the tunneling 
amplitude is modified in a significant way. In the approximation when two independent pairs 
of subbands are included, we obtain an enhanced transmission probability when the barrier 
height is close to TujJq. This is due to the fact that one of the terms in the corresponding wave 
function is perfectly transmitted. By including more than two independent pairs, this effect 
will decrease since perfect transmission will occur only for the two wave function terms. The 
effect is not sensitive to the barrier width, and therefore, can be considered as a reminiscent 
of the Klein paradox. 

We have introduced disorder phenomenologically in this paper through non-conservation 
of the transverse electron momentum component, which is shown to suppress the perfect 
transmission along the diagonal k^i = k^^- From a physical point of view, this disorder 
model could be interpreted as arising from surface roughness of the potential barrier. The 
same type of statistical distribution can be applied to fluctuations in the the barrier width 
which will result in modification of the intensity and location of the transmission peaks. 
Consequently, the transmission maxima observed experimentally will not exactly match to 
those theoretically predicted for clean samples. 
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FIG. 1: (Color online) Transmission probability T of electrons as functions of incoming energy 
e, incident angle (/>, and longitudinal wave numbers hvpkxx and hvpkx- In (a), the transmission 
probability (purple dashed curve) is plotted as a function of the incident particle energy for head- 
on collision {(p = 0). The red solid line is the transmission probability T = 1 (Klein paradox) for 
infinite graphene when the energy gap A = 0. In (b), we show the angular distributions of T for 
infinite graphene (solid red curve) and the dressed states (purple dashed curve) with e = Vo/6. In 
(c) and (d), the transmission probability is plotted as a function of the longitudinal momenta in 
regions "1" and "2" (before the potential barrier and in the barrier region, respectively) for infinite 
graphene (on the left) and irradiated graphene (on the right). The diagonal kxi = kx2 corresponds 
to the absence of a potential barrier (Vq = 0) and displays a perfect transmission with T = 1 
as long as both longitudinal momenta are real determined by both the barrier height Vq and the 

energy gap A. 
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FIG. 2: (Color online) Transmission probability T as functions of both the incoming electron energy 
e and the angle of incidence (p. Plot (a) (far left) is for infinite graphene (A = 0) with obvious 
Klein paradox for = 0. Plots (b) and (c) show the transmission T for electron dressed states, 
in the two nearest subbands approximation, with gap energies A = Vo/15 (middle) and A = Vq/5 
(far right). 
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FIG. 3: (Color online) Energy dispersion and transmission probability for the case of two inde- 
pendent pairs of energy subbands. Panel (a) on the left shows the energy dispersion by taking 
into account the coupling between the states corresponding to two different photon occupation 
numbers in the linear approximation. In the middle panel (b), we present a schematic diagram 
showing how the transmission probability may be increased due to the presence of the other energy 
subbands. Panel (c) on the right gives the transmission vs. incoming particle energy and the angle 
of incidence for the case of two independent pairs of energy subbands. 
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FIG. 4: (Color online) Effect of disorder on the transmission probability using the Lorentzian 
distribution model. For the panel (a) on the left, we plot in (i) the transmission probability as a 
function of k^i for dressed electrons (black curve) compared to infinite graphene (red curve) for 
head-on collision. The plot in (ii) shows a two-peak model distribution under the influence of 
disorder. For panel (b) on the right, we show the transmission probability for electron dressed 
states with A = Vo/3 in the presence of disorder. 
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